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Abstract 

We develop a quantum theory of near-field electro dynamical properties of 
carbon nanotubes and investigate spontaneous decay dynamics of excited 
states and van der Waals attraction of the ground state of an atomic sys- 
tem close to a single-wall nanotube surface. Atomic spontaneous decay ex- 
hibits vacuum-field Rabi oscillations - a principal signature of strong atom- 
vacuum-field coupling. The strongly coupled atomic state is nothing but 
a 'quasi-lD cavity polariton'. Its stability is mainly determined by the atom- 
nanotube van der Waals interaction. Our calculations of the ground-state 
atom van der Waals energy performed within a universal quantum mechani- 
cal approach valid for both weak and strong atom-field coupling demonstrate 
the inapplicability of conventional weak-coupling-based van der Waals inter- 
action models in a close vicinity of the nanotube surface. 

Keywords: Carbon nanotubes, Quantum electrodynamics, Near-field effects 



* To be published in "Nanotubes: New Research", edited by F.Columbus (Nova 

Science, New York, 2005). 
^ Corresponding author. E-mail address: bondarevOtut . by 



Contents 

1 Introductioril 

2 Field Quantization Formalisml 



3 Total Hamiltonian in terms of the local photonic DOS 



3 
6 
10 



4 Soontaneous Decav Dynamics of an Excited Atomic State near a 


Carbon Nanotubd 


4.1 


Qualitative Analysis 


16 




4.1.1 Markoyian AoDroximation 


16 




4.1.2 Sinffle-Resonance ADoroximation of the Local Photonic DQS 


17 


4.2 


Numerical Results and Discussion! 


18 



5 van der Waals Energy of a Ground-State Atom neara Carbon Nanotubd 23 



5.1 Qualitative Analvsi 



26 



5.1.1 Weak Coupling Regime! 26 



5.1.2 Strong Coupling Regime! 27 



5.2 Numerical Results and Discussion! 27 

6 Conclusion! 33 
A Two-level approximation! 36 

' 38 

44 



B Green tensor and local photonic DQS 
C Relation of Eq. ( 70]) w^ith the Casimir Polder formula! 



2 



1 Introduction 



Carbon nanotubes (CNs) are graphene sheets rolled-up into cylinders of approxi- 
mately one nanometer in diameter. Extensive work carried out worldwide in recent 
years has revealed the intriguing physical properties of these novel molecular scale 
wires Nanotubes have been shown to be useful for miniaturized electronic, 

mechanical, electromechanical, chemical and scanning probe devices and materi- 
als for macroscopic composites 0- Important is that their intrinsic properties may 
be substantially modified in a controllable way by doping with extrinsic impurity 
atoms, molecules and compounds 4 . Recent successful experiments on encapsula- 
tion of single atoms into single-wall CNs and their intercalation into single-wall 
CN bundles ||4, ,6| , along with numerous studies of monoatomic gas absorbtion by 
the CN bundles (see [7j for a review), open routes for new challenging nanopho- 
tonics applications of atomically doped CN systems as various sources of coherent 
light emitted by dopant atoms. This, in turn, stimulates an in-depth theoretical 
analysis of near-field electrodynamical properties of atomically doped CNs. Of pri- 
mary importance towards the nanophotonics applications are the problems of atomic 
spontaneous emission and atom-nanotube van der Waals (vdW) interactions. The 
problems are not only of applied but also of fundamental interest as they shed light 
on the peculiarities of the atom-electromagnetic-field interactions in low-dimensional 
dispersive and absorbing surroundings. 

It has long been recognized that the spontaneous emission rate of an excited 
atom is not an immutable property, but that it can be modified by the atomic en- 
vironment. Generally called the Purcell effect [8_, the phenomenon is qualitatively 
explained by the fact that the local environment modifies the strength and distri- 
bution of the vacuum electromagnetic modes with which the atom can interact, 
resulting indirectly in the alteration of atomic spontaneous emission properties. The 
Purcell effect took on special significance recently in view of rapid progress in physics 
of nanostructures. Here, the control of spontaneous emission has been predicted to 
have a lot of useful applications, ranging from the improvement of existing devices 
(lasers, light emitting diodes) to such nontrivial functions as the emission of nonclas- 
sical states of light _9^. In particular, the enhancement of the spontaneous emission 
rate can be the first step towards the realization of a thresholdless laser jJU] or 
a single photon source The possibility to control atomic spontaneous emis- 
sion was shown theoretically for microcavities and microspheres |121 1131 114j , optical 
fibers p^, photonic crystals ll^, semiconductor quantum dots Recent tech- 
nological progress in fabrication of low-dimensional nanostructures has enabled the 
experimental investigation of spontaneous emission of single atoms ^Hl and semi- 
conductor quantum dots jl9l I2UI I21j in microcavities, single molecules in photonic 
crystals |^ , quantum dots in photonic crystal nanocavities |,23i .24] . 

Typically, there may be two qualitatively different regimes of interaction of an 
excited atomic state with a vacuum electromagnetic field in the vicinity of the CN. 



3 



They are the weak couphng regime and the strong coupUng regime |^ . The former 
is characterized by the monotonous exponential decay dynamics of the upper atomic 
state with the decay rate altered compared with the free-space value. The latter is, in 
contrast, strictly non-exponential and is characterized by reversible Rabi oscillations 
where the energy of the initially excited atom is periodically exchanged between the 
atom and the field. A theoretical analysis recently done [261 1271 128j of spontaneous 
decay dynamics of excited atomic states near CNs has brought out fascinating pe- 
culiarities of the vacuum-field interactions in atomically doped CNs. The relative 
density of photonic states (DOS) near a CN effectively increases due to the pres- 
ence of additional surface photonic states coupled with CN electronic quasiparticle 
excitations [22]. This causes an atom- vacuum-field coupling constant, which is pro- 
portional to the photonic DOS, to be very sensitive to the atom-CN-surface distance. 
If the atom is close enough to the CN surface and the atomic transition frequency is 
in the vicinity of a resonance of the photonic DOS, the system shows strong atom- 
vacuum-field coupling giving rise to rearrangement ("dressing") of atomic levels by 
the vacuum-field interaction [271 128j . If the atom moves away from the CN surface, 
the atom-field coupling strength decreases, smoothly approaching the weak coupling 
regime at large atom-surface distances since the role of surface photonic states dimin- 
ishes, causing the relative photonic DOS to decrease. This suggests strictly nonlinear 
atom-field coupling and a primary role of the distance-dependent (local) photonic 
DOS in the vicinity of the CN, so that conventional (weak-coupling) atom-field in- 
teraction models based upon standard vacuum quantum electrodynamics (QED) 
as well as those using the linear response theory (see [291 l3Uj for a review) are, in 
general, inapplicable for an atomic system in a close vicinity of a carbon nanotube. 

The aforesaid is extremely important for a proper understanding the monoatomic 
gas absorbtion processes by the CNs, as they are those arising from the atom- 
nanotube vdW interaction, which, from a fundamental point of view, is the result of 
electrodynamic coupling between the polarization states of the interacting entities 
that originates from the vacuum fiuctuations of an electromagnetic field. Modern 
ab-initio methods based on the density functional theory (DFT) let one study the 
cohesive and adsorption properties of single-wall CNs and CN bundles to an ade- 
quate level of accuracy 31 . For example, recent hydrogen [321 133j and oxygen j^l] 
adsorption DFT-based studies explore the CN gas-storage capacity and the nature 
of hole doping processes of the semiconducting CNs under oxygenation. However, 
DFT is known to be reliable in describing short-range electron correlation effects 
whereas the vdW energy is contributed by both short-range and long-range interac- 
tions p5| • The short-range contribution consists of the repulsive part and attractive 
part coming, respectively, from the overlap of core electrons on adjacent molecules 
and from the decrease in electron kinetic energy due to delocalization. The long- 
range contribution is known as the London energy which originates from the 
(attractive) dispersion interaction of fiuctuating dipoles and dominates the short- 
range contributions at distances larger than the equilibrium one. While adequate 
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in describing the first two contributions, DFT may fail in reproducing the long- 
range dispersion forces correctly, especially in graphitic structures as discussed in 
Ref. |35j . An alternative (classical) approach for computing the London dispersion 
energy is by summing up the empirical pair potentials with parameters fitted 
to experimental data, or for graphitic structures often taken from highly oriented 
pyrolytic graphite [SZl EH] • 

We develop a quantum theory of near-field electro dynamical properties of atom- 
ically doped CNs and apply it to the study of atomic spontaneous emission and 
London-type dispersion atom-nanotube interaction (simply called the vdW inter- 
action throughout this Chapter). Our theory is based on a unified macroscopic 
QED formalism developed for dispersing and absorbing media in Ref. |12| I13j. The 
formalism generalizes the standard macroscopic QED normal-mode technique by 
representing mode expansion in terms of a Green tensor of the (operator) Maxwell 
equations in which material dispersion and absorbtion are automatically included. 
In more detail, the Fourier-images of electric and magnetic fields are considered as 
quantum mechanical observables of corresponding electric and magnetic field op- 
erators. The latter ones satisfy the Fourier-domain operator Maxwell equations 
modified by the presence of a so-called operator noise current written in terms of 
a bosonic field operator and an imaginary dielectric permittivity. This operator is 
responsible for correct commutation relations of the electric and magnetic field op- 
erators in the presence of medium-induced absorbtion. The electric and magnetic 
field operators are then expressed in terms of a continuum set of operator bosonic 
fields by means of the convolution of the operator noise current with the electro- 
magnetic field Green tensor of the system. The bosonic field operators create and 
annihilate single-quantum electromagnetic medium excitations. They are defined by 
their commutation relations and play the role of the fundamental dynamical vari- 
ables in terms of which the Hamiltonian of the composed system "electromagnetic 
field + dissipative medium" is written in a standard secondly quantized form. 

The Chapter is organized as follows. In Section [3 using the approach sum- 
marized above, we describe the general electromagnetic field quantization scheme 
and derive the total Hamiltonian for an atomic system (an atom or a molecule) 
nearby a single-wall carbon nanotube. The nanotube is considered as an infinitely 
long, infinitely thin, anisotropically conducting cylinder. Its surface conductivity is 
represented in terms of the 7r-electron dispersion law obtained in the tight-binding 
approximation with allowance made for the azimuthal electron momentum quanti- 
zation and the axial electron momentum relaxation (HHI • Only the axial conductivity 
is taken into account and the azimuthal one, being strongly suppressed by trans- 
verse depolarization fields |^ |^ HH UHl ll^i is neglected. In Section |31 based 
on the total Hamiltonian obtained in Section [2 we derive the simplified secondly 
quantized Hamiltonian representing the "atom-nanotube" coupled system in terms 
of two standard approximations. They are the electric dipole approximation and 
the atomic two-level approximation. This Hamiltonian generalizes that earlier ap- 
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plied for the atomic spontaneous decay analysis near spherical microcavities |12| 
and nanotubes j27l I28j in that the rotating wave approximation is not used and the 
diamagnetic term of the atom-field interaction is not neglected. This Hamiltonian is 
a starting point to derive an evolution equation for the population probability of the 
upper state (Section and a van der Waals energy equation of the lower (ground) 
state of the two-level atomic system (Section [SJ coupled with the CN modified 
vacuum electromagnetic field. The equations are represented in terms of the local 
photonic DOS and are valid for both strong and weak atom-field coupling. In Sec- 
tion [IJ we demonstrate that when the atom is close enough to the nanotube surface 
and the atomic transition frequency is in the vicinity of a resonance of the local 
photonic DOS, the upper state of the coupled system decays via underdamped Rabi 
oscillations which are a principal signature of strong atom-vacuum-field coupling 
- the result recently detected for quasi-two-dimensional (2D) excitonic [IHI, inter- 
subband electronic [1^1 and quasi-OD excitonic jl9l I21j transitions in semiconductor 
quantum microcavities as well as for quasi-OD excitonic transitions in photonic crys- 
tal nanocavities |23| . In Section |SJ to take strong atom- field coupling into account in 
a correct way and thereby to advance in physical understanding the atom-nanotube 
vdW interactions deeper than simple analyzing empirical model potentials j3S]) we 
develop a universal quantum mechanical approach to the atom-nanotube vdW en- 
ergy calculation. The approach was first suggested in Refs. [13 EE]- It is based on 
the perturbation theory for degenerate atomic levels (see, e.g., [HI), thereby allow- 
ing one to account for both strong and weak atom-vacuum-field coupling regimes in 
a simplest way. Within this approach, the ground-state atom vdW energy is given 
by the integral equation which reproduces a well-known perturbation theory result 
in the weak coupling regime. By solving it numerically, we demonstrate the inap- 
plicability of conventional weak-coupling-based vdW interaction models in a close 
vicinity of the nanotube surface. A summary and conclusions are given in Section|EJ 



2 Field Quantization Formalism 

The quantum theory of the near-field electro dynamical properties of carbon nan- 
otubes involves an electromagnetic field quantization procedure in the presence of 
dispersing and absorbing bodies. Such a procedure faces difficulties similar to those 
in quantum optics of 3D Kramers-Kronig dielectric media where the canonical quan- 
tization scheme commonly used does not work since, because of absorption, corre- 
sponding operator Maxwell equations become non-Hermitian [SIl- As a consequence, 
their solutions cannot be expanded in power orthogonal modes and the concept of 
modes itself becomes more subtle. We, therefore, use a unified macroscopic QED 
approach developed in Refs. jl2lll3j . In this approach, the Fourier-images of electric 
and magnetic fields are considered as quantum mechanical observables of corre- 
sponding electric and magnetic field operators. The latter ones satisfy the Fourier- 
domain operator Maxwell equations modified by the presence of a so-called operator 
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Figure 1: The geometry of the problem. 



noise current density J(r, cu) written in terms of a 3D vector bosonic field operator 
f(r,a;) and a medium dielectric tensor e(r,a;) (supposed to be diagonal) as^ 

Jj(r,u;) = ■^y^hlmeii{r,uj) fi{r,uj), i = 1,2,3. (1) 
zvr 

This operator is responsible for correct commutation relations of the electric and 
magnetic field operators in the presence of medium-induced absorbtion. In this 
formalism, the electric and magnetic field operators are expressed in terms of a 
continuum set of the 3D vector bosonic fields f(r,a;) by means of the convolution 
over r of the current with the classical electromagnetic field Green tensor of 
the system. The bosonic field operators F(r,a;) and f(r,a;) create and annihilate 
single-quantum electromagnetic medium excitations. They are defined by their com- 
mutation relations 

[/,(r,a;),/j(r',a;')] = 5,,5{r - r')5{io - co') , (2) 
[/,(r,^),/,(r',u;')] = [ //(r, ^), /j(r', ] = 

and play the role of the fundamental dynamical variables in terms of which the 
Hamiltonian of the composed system "electromagnetic field + dissipative medium" 
is written in a standard secondly quantized form as 

H = J dioniuJdri\r,uj)-i{r,uj). (3) 

Consider a neutral atomic system (an atom or a molecule) with its centre-of-mass 
positioned at the point near an infinitely long single-wall CN. Since the problem 
has a cylindric symmetry, assign the orthonormal cylindric basis {er,e^,ez} with 

^We use Gaussian units. Conversion tables to SI can be found elsewhere (see, e.g., |51p. 
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directed along the nanotube axis and = r^e^ = {rA,0,0} (see Fig. P). For 
carbon nanotubes, their strong transverse depolarization along with transverse elec- 
tron momentum quantization allow one to neglect the azimuthal current and radial 
polarizability ^39i i40l I^Tl 1421 UHl I44j . in which case the dielectric tensor components 
err and e^^ are identically equal to unit. The component e^^ is caused by the CN 
longitudinal polarizability and is responsible for the axial surface current parallel 
to the vector. This current may be represented in terms of the ID bosonic field 
operators by analogy with Eq. Q. Indeed, taking into account the dimensionality 
conservation in passing from bulk to a monolayer in Eq. (jSJ and using a simple 
Drude relation j41j 

aUK,u;) = -zu; ^^^^^ , (4) 

where R ={Rcn, (p, Z} is the radius-vector of an arbitrary point of the CN surface, 
Ren is the radius of the CN, (T^2(R, w) is the CN surface axial conductivity per unit 
length, S is the area of a single nanotube, px is the cubic density of the tubule, one 
immediately has from Eq. 



J(R,a;) = Y f{Ti,uj)ez (5) 

with /(R, w) being the ID bosonic field operator defined on the CN surface. The 
operators P(R,uj) and /(R, w) satisfy the commutation relations on the surface 
of the CN. They create and annihilate the single-quantum bosonic-type electromag- 
netic medium excitation of the frequency lo at the point R of the CN surface. 

We are now in a position to define the Hamiltonian of the system. Following 
the general macroscopic QED approach in the presence of absorbing and dispersive 
media jl3j , we introduce the total nonrelativistic Hamiltonian of an atomic system 
interacting with a CN modified quantized vacuum electromagnetic field in the form 

oo 



H= I dujhw j d'Rf\n,uj)f{n,u) (6) 

Pi - ^A(rA + f-j) ^\ j dr pA{r)(pA{r) + j dv pA{r)Lp{Y) , 



1 

2m- 



where mj, qi, fj and pj are, respectively, the masses, charges, coordinates (relative to 
r^) and momenta of the particles constituting the atomic subsystem. The first term 
of the Hamiltonian describes the quantized medium- assisted vacuum electromagnetic 
field in the presence of the CN. The second and the third terms represent the kinetic 
energy of the charged particles and their mutual Coulomb interaction, respectively, 
with 

V?A(r) = / dv (7) 
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being the scalar potential of the charged particles distributed with the density 
/5/i(r) = qi5{r — r^i — f j) in the atomic subsystem. The last term of the Hamilto- 
nian accounts for the Coulomb interaction of the particles with the CN. The vector 
potential A and the scalar potential of the CN modified electromagnetic field are 
given for an arbitrary r = {r, ip, z} in the Schrodinger picture by 

r°° c -± 

A(r) = dLO—B (r,u;)+h.c., (8) 
Jo 

-V(^(r) = / dwE"(r,u;) +h.c., (9) 
Jo 



where 



E (r, 



^cj) = y"dr'<5^(ll)(r-r') •E(r',w) (10) 
is the transverse (longitudinal) electric field operator with 

6l,{r) = -V.V,^ (11) 

and 

6^f,{r)=6a^6{r)-6lpir) (12) 

being the longitudinal and transverse dyadic (5-functions, respectively, and E repre- 
senting the total electric field operator which satisfies the following set of Fourier- 
domain Maxwell equations 

V x'E{r,uj) =ikU{r,uj), (13) 



V x'H{r,uj) = -ikE{r,Lo) + — I{r,Lo). (14) 

c 

Here, H stands for the magnetic field operator, k = lo/c, and 

I(r,u;)= y"dR5(r-R)J(R,w) =2J(i?,„,(/j,z,w)<5(r-i?en), (15) 

is the exterior operator current density [with J(R, cj) defined by Eq. ©] associated 
with the presence of the CN. 

Prom Eqs. (O-lITni) in view of Eq. ©, it follows that 

E(r,cj) =i—k /(iRG(r,R,a;)-J(R,cj) (16) 

[and H = (iA;)^^V x E accordingly], where G is the Green tensor of the classical 
electromagnetic field in the vicinity of the CN. The Green tensor components satisfy 
the equation 

^ (V X V X - k\^ Gazir, R, c^) = ,5(r - R), (17) 

a=r,ip,z 
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together with the radiation conditions at infinity and the boundary conditions on 
the CN surface. The set of Eqs. (|^- (|17() forms a closed electromagnetic field quan- 
tization formalism in the presence of dispersing and absorbing media which meets 
all the basic requirements of the standard quantum electrodynamics jl3| . All the 
information about medium, which is a carbon nanotube in our case, is contained 
in the Green tensor G via the CN surface axial conductivity Ozz in Eq. which 
comes into play when G is imposed the boundary conditions on the CN surface. 
The classical electromagnetic field Green tensor of the " atom-nanotube" system is 
derived in Appendix B. 

Assuming further that the atomic subsystem is sufficiently localized in space, so 
that the long-wavelength approximation applies, one can expand the field operators 
A(r) and <^(r) in the Hamiltonian © around the atomic center of mass position 
and only keep the leading non- vanishing terms of the expansions. Then, under the 
condition of the Coulomb gauge [pi,A] = 0, one arrives at the approximate total 
Hamiltonian of the form 

H = Hf + Ha + H^aI + H^AF ' (18) 

where 

dujhujJdRp(R,uj)f(R,uj), (19) 



^-^ 2mi 

I Kj 



H^F = - E — P^- A(rA) + d-V^ivA) , (21) 



rriiC 

2m jC 



are, respectively, the Hamiltonian of the vacuum electromagnetic field modified by 
the presence of the CN, the Hamiltonian of the atomic subsystem, and the electric 
dipole approximation for the Hamiltonian of the atom-field interaction (separated 
into two contributions according to their role in the vdW energy - see below) with 
d = J2i Qi^i being the electric dipole moment operator of the atomic subsystem. 



3 Total Hamiltonian in terms of the local photonic DOS 

Starting from Eqs. (HHJ-lESl) and using Eqs. (jHll-jni), and one obtains 
the following two-level approximation for the total Hamiltonian of the system (see 
Appendix A) 

H = Ho + Hint , (23) 
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where the first term stands for the unperturbed Hamiltonian given by 

Ho = Hp + Ha (24) 
with Hp being the Hamiltonian (|19|) of the CN modified field and 

nA = ^a, (25) 

representing the 'effective' unperturbed two-level atomic subsystem, and the second 
term stands for their interaction given by 

dujjdRl g(+) (fa, R, w) - g(-) (rA, R, w) a ] /(R, cu) + h.c. (26) 

Here, the Pauli operators 

a = \l){u\, at = |n)(/|, a, = \u) {u\ - \l) {l\ , (27) 
\u){u\ + \l){l\ = i, 

describe electric dipole transitions between the two atomic states, upper \u) and 
lower |/), separated by the transition frequency iva- This 'bare' transition frequency 
is modified by the interaction (|22|) which, being independent of the atomic dipole 
moment, does not contribute to mixing the |m) and |/) states, giving rise, however, 
to the new renormalized transition frequency 
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{hwA)'' 



/"OO /" 

j du j (iR|g"^(r^,R, u;) 



(28) 



in the 'effective' atomic Hamiltonian (|25j) . On the contrary, the interaction (|21j) . 
being dipole moment dependent, mixes the \u) and |/) states, yielding the pertur- 
bation Hamiltonian (|26|) with the (dipole) interaction matrix elements of the form 

g(±)(rA,R,a;) =g^(rA,R,u;)±— gll(rA,R,t^), (29) 

LOA 

where 

g^^\\\vA,n,uj) = -i^ d.y/nhw Re a,,(R, uj) ^("^^..(rA, R, w) (30) 

with dz = {l\dz\u) = {u\dz\l) being the matrix element of the atomic dipole moment 
z-component and 

^(ll)G,,(rA,R,a;) = J dr6jP{rA-r)Gzzir,n,iv) (31) 

representing the field Green tensor zz-component transverse (longitudinal) with 
respect to the first variable. This is the only Green tensor component we have to 
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take account of. All the other components can be safely neglected, because our 
model neglects the azimuthal current and radial polarizability of the CN from the 
very outset [see the discussion above Eq. (jlj)]. 

Matrix elements H29|) and ()3U() have the following properties (Appendix B) 

J dK Ig^(ll) (r^, R, = ro(u;) ^^^D (r^, c) (32) 



and 

(33) 



y"dR|gW(r^,R,.;)P = i^ro(^) 



Here, .^"'"(ryi, w) and ^"(ryi,^;) are the transverse and longitudinal local (position- 
dependent) photonic DOS functions defined by 

(34) 

ImG^2(r^,r^,a;) 

where 

^(ll)Gi-ill)(rA,rA,u;) = J dr dr' 6i^\\\r a - r) G,,(r, r', u;) .J^dD (r' - r^) (35) 

is the Green tensor zz-component transverse (longitudinal) with respect to both 
coordinate variables, 

ImG°^(rA,rA,t^) = - — (36) 
one 



is the vacuum imaginary Green tensor zz-component j52j. and 

ro(^) = ^^ImGL(rA,rA,^) = (37) 

is the atomic spontaneous decay rate for the z-oriented atomic dipole in vacuum, 
taken at an arbitrary frequency oj [SHllMj . 

Note that the transverse and longitudinal imaginary Green tensor components 
in Eq. can be written as 

Im^Gf, = Im Gl, + Im^G^ , (38) 
Im"G!|^ = Im"G]'j, (39) 

with -'-dl^G^J"^ representing the "pure" CN contribution to the total imaginary Green 
tensor. The longitudinal imaginary Green tensor in Eq. (|39|) is totally contributed 
by a (longitudinal) CN static polarization field and, therefore, does not contain the 
vacuum term which is the transverse one by its definition as there are no polarization 
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Coulomb sources in vacuum [^E^J. Thus, Eq. (|5Hl may, in view of Eqs. (|55|l and 
1)391) . be rewritten in the form 



e^(r^,a;) = 1 + ^ (rA,c^), (40) 
C"(rA,w) = c\rA,uj), (41) 

where the r^-dependent terms come from the presence of the CN. For va > Rem 
they are exphcitly given by (Appendix B) 

j±, ^ 3^ f dys{Raa,x)v{yYll[v{y)u{RMKl[v{y)u{rA)x] 



TT 

p=— oo' 



7h 1 = ^T V f dy siRcn,x)yMy?IpHyMR^n)x]K^[v{y)u{rA)x] 



TT 

P= 



where 



x=^u; (44) 
270 

is the dimensionless frequency with 70 = 2.7 eV being the carbon nearest neighbor 
hopping integral ^3 appearing in the CN surface axial conductivity in Eq. ©, 
Ip and Kp are the modified cylindric Bessel functions, v{y) = y'^ — 1 , u{r) = 
27or/fic, and s{Rcn,x) = 2iau{Rcn)xazziRcm x) with a^z = 27rhazz/e'^ being the 
dimensionless CN surface conductivity per unit length and a = e^/hc = 1/137 
representing the fine-structure constant. The integration contour C runs along the 
real axis of the complex plane and envelopes the branch points y = ±1 of the function 
v{y) in the integrands from below and from above, respectively. For < Ren, 
Eqs. (|^^ and are modified by a simple replacement <^ Ren in the Bessel 
function arguments in the numerators of the integrands. Note the divergence of the 
local photonic DOS functions (,'^^^^\rA, x) at va = Ren, i- e. when the atom is 
located right on the CN surface. Mathematically, this comes from the logarithmic 
divergence of the Bessel function Kq at small arguments (see, e.g., ISGj). yielding the 
divergent product IqKq and, as a consequence, the divergent integrands in Eqs. (|42)) 
and ()43() , if and only if the arguments of the functions Iq and Kq in the numerator of 
the integrand vanish simultaneously. Physically, the point is that the CN dielectric 
tensor longitudinal component €zz [which, according to Eqs. (@J), © and (|TH|) is 
responsible for the surface axial conductivity azz in Eqs. (|^^ and is obtained 
as a result of a standard procedure of physical averaging a local electromagnetic field 
over the two spatial directions in the graphene plane [SJ- Such averaging does not 
assume extrinsic atoms on the graphene surface. To take them into consideration 
the averaging procedure must be modified. Thus, the applicability domain of our 
model is restricted by the condition 

\rA- Ren\> a, (45) 
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where a = 1.42 A is the graphene interatomic distance |55j . 

Eqs. H23|) - H3U() represent the total secondly quantized Hamiltonian of an atomic 
system interacting with the electromagnetic field in the vicinity of the carbon nan- 
otube. In deriving this Hamiltonian, there were only two standard approximations 
done. They are the electric dipole approximation and the two-level approximation. 
The rotating wave approximation commonly used was not applied, and the (dia- 
magnetic) A^-term of the atom-field interaction [Eq. ((221)] '^^s not neglected. The 
latter two approximations are justified in analyzing an atomic state evolution due 
to real dipole transitions to the nearest states (shown schematically in Fig. |21) |2S| , 
described for two-level systems by the g'^"*")-terms of the Hamiltonian H26|) . They 
are, however, inappropriate in calculating the energy, where the initial energy level 
is being shifted due to the perturbation caused by virtual dipole transitions via 
intermediate atomic states (see Fig. El in Sectional), so that the system always re- 
mains in the initial state and the terms neglected contribute to the level shift to the 
same order of magnitude. For two-level systems, this is described by g(~)-terms of 
the Hamiltonian ()26() (usually neglected within the rotating wave approximation) 
together with the atomic transition frequency H28|l renormalized due to the (usually 
discarded) A^-term of the atom-field interaction. 



4 Spontaneous Decay Dynamics of an Excited Atomic 
State near a Carbon Nanotube 

When the atom is initially in the upper state and the field subsystem is in vacuum 
(see Fig. 121), the time-dependent wave function of the whole system can be written as 

\m) = Cuit)e-'^'^^/''^'\u)m) 

/roo 
dKJ dioCiiR,LO,t)e-'^'^-^^/^^^\l)\{l{K,oj)}), (46) 

where |{0}) is the vacuum state of the field subsystem, |{1(R, u;)}) is its excited state 
where the field is in a single-quantum Fock state, Cu and Ci are the population 
probability amplitudes of the upper state and lower state of the whole system, 
respectively. The function (|46)) accounts for the lower state degeneracy with respect 
to R as well as its possible degeneracy with respect to u! when to ~ uja as is seen in 
Fig. 121 

The time-dependent Schrodinger equation with the Hamiltonian (|23j) - (|2H)) and 
the wave function ()46() yields the following evolution law for the population proba- 
bility amplitude Cu of the upper state 

Cu{t) = 1 + /V(t - t') Cu{t') dt', (47) 
Jo 

K{t-t') = T2 du^— — /dR|g(+)(rA,R,a;)p. (48) 
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Figure 2: Schematic of the energy levels of the "atom-nanotube" coupled system in 
the problem of the spontaneous decay. On the left are the unperturbed atomic levels 
given by the Hamiltonian (|25|) . on the right are the levels of the coupled system. 
The relevant real transition of the system is shown by the vertical dashed arrow. 
See Eqs. ^ and (gH) for uja- 



This, upon recalling Eqs. (1331) , (j^H)) . (|lT|l and (|1H) and introducing the dimensionless 
time 

270 



reduces in dimensionless variables to 

Cu{t) = 1 + rK{T - t') Cu{r') dr', 



(49) 



(50) 



K{t - t') = ^ 



^2 roo 



dx 



fo(x) 



1 + e {rA,x) + 



-)'f{rA,x) 
xaJ 



-i{x-XA){T-T') _ 

i {x - xa) 



(51) 



with the r^-dependent transverse and longitudinal local photonic DOS functions 
(,'^{rA,x) and ^"(r^,a;) given explicitly by Eqs. (|^^ and respectively, 



h 



Adl /27oV 



X 



ro(x) = — ro(x)-^^^3 ^ ^ J 

being the dimensionless vacuum spontaneous decay rate H37() . and 



(52) 



XA = Xa 



vr 



dx 



(53) 



representing the dimensionless renormalized transition frequency H28|l where a di- 
vergent contribution to the vacuum Lamb shift has been omitted since the true 
(renormalized) vacuum Lamb shift may be considered as included in the atomic 
transition frequency xa ^11 • 
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Eq. H5U|) is a well-known Volterra integral equation of the second kind, which 
in our case describes spontaneous decay dynamics of an excited two-level atomic 
system interacting with the vacuum electromagnetic field modified by the presence 
of a CN. All the CN parameters that are relevant for the spontaneous decay are 
contained in the local photonic DOS functions in the kernel H51|). The latter ones 
are, in view of Eqs. H34|) . ()4U|) and (|41|) . determined by the imaginary classical Green 
tensor of the CN modified electromagnetic field. 



4.1 Qualitative Analysis 

Let us first analyze qualitatively the time dynamics of the upper state population 
probability Cu{t) in terms of two different approximations admitting analytical 
solutions of the evolution problem (|50|) , (|51j) . They are the Markovian approximation 
and the single-resonance approximation of the local photonic DOS in the vicinity of 
the CN. 



4.1.1 Markovian Approximation 



Suppose the atom is faraway from the CN. Obviously, the local photonic DOS 
functions ^"""(r^jx) and ^"(r^,x) are then small in their values and the atom-field 
coupling matrix elements H32|) and (|33() are mainly contributed by the free space 
term which comes from the unity in [see Eqs. (jSU, (|!-{8j) ~ (|41jl ]. In this particular 
case, the Markovian approximation is applicable for sure. The latter one is known to 
be an excellent approximation for describing the radiative decay of an excited atom 
in free space where the atom-field coupling strength is weak enough for atomic 
motion memory effects to be insignificant, so that they may be neglected [HSI- The 
time-dependent factor in the kernel (|51() may then be replaced by its long-time 
approximation (see, e.g., [I9] ) 

— > -t:5{x - xa) + i — ■, (54) 

i{x — xa) X — xa 

{V denotes a principal value) where xa ~ xa according to Eq. with small enough 

^ (r^,x). In so doing, the kernel becomes 

1 



K{t -t') = - 
and the evolution equation (|5Ujl yields 



■r(xA) +iA{xA) 



Cu{t) = exp 



--f(xA) +iA(xA) 



the exponential decay dynamics of the upper atomic level shifted by 

2 



ax 



XA 



^) t{rA,x)+i\vA,x) 



(55) 



(56) 



(57) 
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[here, a divergent contribution to the vacuum Lamb shift has been omitted for the 
same reason as that in Eq. (|53|) ] with the rate 



t{xA)=ro{xA) [l+^^irA,XA)+^hrA,XA)] (58) 
accounting for the presence of the nanotube. 

4.1.2 Single-Resonance Approximation of the Local Photonic DOS 

Another approximation that admits an analytical solution of the evolution problem 
(jSnj, H51() is the single-resonance approximation of the local photonic DOS. Suppose 

the functions (,^{rA,x) and ^"(r^,x) have sharp peaks at x = Xr [the fact that 
they have the same positions of extrema is obvious from Eqs. (|42|) and (|4l-ij) : see 
also the footnote on page I28j . For all x in the vicinity of Xr, their shape may then 
be roughly approximated by the Lorentzian functions with the same half-width-at- 
half-maximum 6xr of the form 

e {rA,x)^c (r^'^^)(^_^^)2+5^2- m 

The integrand in the kernel (|51|) is seen to be basically contributed hy x ^ xa, 
thus making it possible to approximate the factor in the square brackets as follows 

l+l^irA,x) + (^)'f{rA,x) « 1 + C^(r^,x) + (^)'f{rA,x) 
\XAj \XAj 

with XA given in an explicit form by Eq. (|53() . Keeping this and the Lorentzian 
approximation ()59() in mind, one calculates the kernel (|51l) analytically to linear 
approximation in 5xr to obtain 

K{r - r') ^ f (X.) ^ ^"P[-^("'- - '^^r- ~XA){r- r')] - 1 ^ ^ ^, ^^^^ 
2 iyxr — ioxr — xa) 

with V{xr) given by Eq. (|5H)) . Substituting this into Eq. (fSU)) and making the differ- 
entiation of both sides of the resulting equation over time, followed by the change 
of the integration order and one more time differentiation, one straightforwardly 
arrives at a second order ordinary homogeneous differential equation of the form 

Cu{r)+i{xr-i5xr-iA)Cu{r) + {X/2fCu{T) = Q, (61) 

where 



X = \j25XrT{Xr) , (62) 
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with the solution given for xa ~ Xr by 



6xr — \/ 5x1 — X'^j — 



5xr + yjdxl - xA ^ 



(63) 



This solution is approximately valid for those atomic transition frequencies xa which 
are located in the vicinity of the photonic DOS resonances whatever the atom-field 
coupling strength is. In particular, if (X/Jxr.)^ <C 1, Eq. (|63|) yields the exponential 
decay of the upper atomic state population probability |Cu(r)p with the rate T{xr) 
in full agreement with Eq. ()56() obtained within the Markovian approximation for 
weak atom-field coupling. In the opposite case, when {X/5xr)'^^l, one has 



\Cu{r) 



e-^^^"cos2 



'Xt 



(64) 



and the decay of the upper atomic state population probability proceeds via damped 
Rabi oscillations with the dimensionless Rabi frequency X given by Eq. (|62)) . This 
is the principal signature of strong atom-field coupling which is beyond the Marko- 
vian approximation. Expressing T{xr) in Eq. (|62|) in terms of the local photonic DOS 
functions by means of Eqs. (jSHj) . (|4U() . H41|) and approximating for atomic systems 
with Coulomb interaction Tq{x) ~ 137~^x |49j . one may conveniently rewrite the 
strong atom-field coupling condition in the form 



2x^ 



i^{Xr)+iKxr) 
1373 5Xr 



» 1, 



(65) 



from which it follows that the strong coupling regime is fostered by high and narrow 
resonances in the local photonic DOS. 



4.2 Numerical Results and Discussion 

To get beyond the Markovian and single-resonance approximations, we have solved 
Eqs. (fSO]) and numerically. The exact time evolution of the upper state popula- 
tion probability |Cu(t)P was obtained for the atom placed [in a way that Eq. (|45|) 
was always satisfied] in the center and near the wall inside, and at different distances 
outside achiral CNs of different radii. The local photonic DOS functions were com- 
puted according to Eqs. (|in|) - H43|) . The CN surface axial conductivity Gzz appearing 
in Eqs. (|42j) . (|43j) was calculated in the relaxation-time approximation with the relax- 
ation time 3 x 10"^^ s; the spatial dispersion of azz was neglected 3!) . The vacuum 
contribution to the kernel (|51|) (that coming from the unity in the square brackets) 
was taken into account within the Markovian approximation (|54() with the divergent 
vacuum Lamb shift omitted as the true (renormalized) vacuum Lamb shift may be 
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considered as included in the 'bare' atomic transition frequency xa jl2l I13j . The 
Markovian approximation is known to be an excellent approximation for describing 
atom-radiation interaction processes in free space j58| . The vacuum spontaneous 
decay rate was approximated by the expression FoC^c) ~ 137~^x valid for atomic 
systems with Coulomb interaction j49j . 

Figure 01 (a) presents ^"""(x) for the atom in the center of the (5,5), (10,10) and 
(23,0) CNs. It is seen to decrease with increasing the CN radius, representing the 
decrease of the atom-field coupling strength as the atom moves away from the CN 
surface. To calculate |C„(t)P in this particular case, we have fixed xa = 0.45 (indi- 
cated by the vertical dashed line), firstly, because this transition is located within 
the visible light range 0.305 <x < 0.574, secondly, because this is the approximate 
peak position of ^"""(x) for all the three CNs. The functions |Cu(t)P calculated are 
shown in Figure 121 (b) in comparison with those obtained in the Markovian approx- 
imation yielding the exponential decay. The actual spontaneous decay dynamics is 
clearly seen to be non-exponential. For the small radius (5,5) CN, Rabi oscillations 
are observed, indicating a strong atom-field coupling regime related to strong non- 
Markovian memory effects. Eq. (|65() is satisfied in this case. With increasing the 
CN radius, as the value of ^^(0.45) decreases, Eq. (|65|) ceases to be valid and the 
decay dynamics approaches the exponential one with the decay rate enhanced by 
several orders of magnitude compared with that in free space. Note that, though 
the distance from the atom to the CN surface is larger for the (23,0) CN than for 
the (10,10) CN, the deviation of the actual decay dynamics from the exponential 
one is larger for the (23,0) CN. This is an obvious consequence of the influence of 
a small neighboring peak in the (23,0) CN photonic density of states [Figure |31 (a)]. 

In Ref. 0, formation of Cs-encapsulating single- wall CNs was reported. In a par- 
ticular case of the (10,10) CN, the stable Cs atom/ion position was observed to be at 
distance of 3 A from the wall. We have simulated the spontaneous decay dynamics 
for a number of typical atomic transition frequencies for this case. Figure01(a) shows 
the local photonic DOS and the five specific transition frequencies xa (dashed lines) 
for which the functions |Cu(t)p presented in Figure |ll(b) were calculated. Rabi os- 
cillations are clearly seen to occur in the vicinity of the highest peak (xa ~ 0.22) 
of the photonic density of states. Important is that they persist for large enough 
detuning values Xyi « 0.21-^0.25. For xa = 0.30, the density of photonic states has 
a dip, and the decay dynamics exhibits no Rabi oscillations, being strongly non- 
exponential nevertheless. For xa = 0.45, the intensity of the peak of the photonic 
density of states is not large enough and the peak is too broad, so that the strong 
atom-field coupling condition (|65|l is not satisfied and the decay dynamics is close 
to the exponential one. 

Figure El (a) shows the local photonic DOS for the atom outside the (9,0) CN 
at the different distances from its surface. The vertical dashed lines indicate the 
atomic transitions for which the functions |Cu(r)p in Figures [51 (b), (c), and (d) 
were calculated. The transitions x^ = 0.33 and 0.58 are the positions of sharp peaks 
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X 



Figure 3: Transverse local photonic DOS's (a) and upper-level spontaneous decay 
dynamics (b) for the atom in the center of different CNs. The ('bare') atomic tran- 
sition frequency is indicated by the dashed line in Fig. |31 (a). 
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Figure 4: (a) Fragment of the transverse local photonic DOS for the atom inside the 
(10,10) CN at distance of 3 A from the wah (the situation observed experimentahy 
for Cs in Ref. 5 ). (b) Upper-level spontaneous decay dynamics for different ('bare') 
atomic transition frequencies [indicated by the dashed lines in Fig. (a)] in this 
particular case. 
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Figure 5: (a) Transverse local photonic DOS's for the atom located at different 
distances outside the (9,0) CN. (b, c, d) Upper-level spontaneous decay dynamics 
for the three ('bare') atomic transition frequencies [indicated by the dashed lines in 
Fig. 13 (a)] at different atom-nanotube-surface distances. 
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(at least for the shortest atom-surface distance), while xa = 0.52 is the position of 
a dip of the function ^-^{x). Very clear underdamped Rabi oscillations are seen for 
the shortest atom-surface distance at xa = 0.33 [Figure 0] (b)], indicating strong 
atom-field coupling with strong non-Markovity. For x^ = 0.58 [Figure |21 (d)], the 
value of ^-"-(O.SS) is not large enough for strong atom-field coupling to occur, so 
that Eq. (|65() is not fulfilled. As a consequence, the decay dynamics, being strongly 
non-exponential in general, exhibits no clear Rabi oscillations. For xa = 0.52 [Fig- 
urelSl(c)], though ^-'-(0.52) is comparatively small, the spontaneous decay dynamics 
is still non-exponential, approaching the exponential one only when the atom is far 
enough from the CN surface. 

The reason for non-exponential spontaneous decay dynamics in all the cases con- 
sidered is similar to that taking place in photonic crystals |16) . When the atom is 
close enough to the CN surface, an absolute value of the relative density of photonic 
states is large and its frequency variation in the neighborhood of a specific atomic 
transition frequency essentially infiuences the time behavior of the kernel (|51|) of 
the evolution equation H50|). Physically, this means that the correlation time of the 
electromagnetic vacuum is not negligible on the time scale of the evolution of the 
atomic system, so that atomic motion memory effects are important and the Marko- 
vian approximation in the kernel H51|) is inapplicable. The drastic increase of the 
local photonic DOS near a CN can be interpreted as being due to the electromag- 
netic vacuum renormalization: the relative density of photonic states near the CN 
effectively increases since, along with ordinary free photons, there appear surface 
photonic states coupled with CN electronic quasiparticle excitations. These latter 
ones are responsible for the nonradiative atomic decay (photon emission by the atom 
with subsequent CN quasiparticle excitation) 26 . Obviously, such a decay does not 
respond sensitively to the microscopic radiation-field structure in a close vicinity 
of the CN surface. This, along with Eq. H45|) . justifies rather small atom-surface 
distances we have considered within the macroscopic Green-tensor approach. 



5 van der Waals Energy of a Ground- State Atom near 
a Carbon Nanotube 

When the atom is in the ground state, its vdW energy near the nanotube is given 
by the r^i-dependent contribution to the ground-state eigenvalue of the total Hamil- 
tonian p3j) - (|26)) . Important, however, is that the unperturbed atomic subsystem 
is now described by the Hamiltonian (|25() with the renormalized transition fre- 
quency (|28j) . Eqs. (|53j) and (|44|1 represent the latter one in terms of the transverse 
r^-dependent photonic DOS ^ (r^, x). They indicate that the transition frequency 
UA decreases with increasing ^ , i. e. when the atom approaches the CN surface 
[see Fig. |21 (a) as an example] , thereby bringing the unperturbed atomic levels to- 
gether, or even making them degenerated if ^ is large enough [see the schematic 
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Figure 6: Schematic of the energy levels in the problem of the ground-state atom 
vdW interaction with the CN. On the left are the unperturbed atomic levels given 
by the Hamiltonian (|25() . on the right are the levels of the coupled "atom-nanotube" 
system given by the total Hamiltonian H23 |) - H26() . The relevant virtual transition of 
the system is shown by the vertical dashed arrows. See Eqs. ()53() and (|44|) for uja- 



of the energy levels of the problem in Fig. El and a typical example of the (9,0) 
CN in Fig. [7j. If one uses the standard perturbation theory, this would yield the 
divergence of high-order corrections to the vdW energy at small atom-CN-surface 
distances because of the small energy denominators in the corresponding power- 
series expansions. To account for a possible degeneracy of the unperturbed atomic 
levels in a correct way, one has to use the perturbation theory for degenerated atomic 
levels (see, e.g., UHl)- In so doing, one has to take the upper state degeneracy into 
account of the whole system with respect to R as well as its possible degeneracy 
with respect to w when lo ~ uja ~ as is seen in Fig. El Thus, the ground-state 
wave function of the whole system should be represented as a coherent mixture of 
the lower and upper states of the form 

/■oo f 

IV) = Q|/)|{0}) + J^dcoJdKCu(R,u;)\u)\{l(R,u;)}), (66) 

where Ci and Cu(R, w) are unknown mixing coefficients of the lower and upper 
states of the whole system. 

The total energy E of the ground state is now given by the solution of a sec- 
ular equation obtained by applying the Hamiltonian H23 |) - H26|) to the wave func- 
tion (|66|) . This yields the integral equation 

2 

which the ground-state-atom vdW energy £'^tf,(ryi) is determined from by means of 
an obvious relationship 

E = -^ + E,^{rA) (68) 
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Figure 7: Dimensionless unperturbed energy levels [given by the Hamiltonian (|24|) 
and (Uni)] and total ground-state energy [given by Eqs. (jSHIl-lIZni), (EH) and (HSJ] as 
functions of the atomic position for the two-level atom outside the (9,0) CN. See 
Eq. (031) for XA- 



with the first term representing the 'bare' (non-interacting) two-level atom in the 
ground-state. This, upon recalling Eqs. (jlU)) . (|1T|) and (|^H) and introducing 

the dimensionless vdW energy 



£vw{t^a) 



Evw{rA) 
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reduces to the integral equation of the form 



XA 
2lT 



dx 



(69) 



(70) 



^1 



fo(x) 



^ {rA,x) + 



\2 



C (ryi,x) 



re -I- XA 



27r 



°^ ro(x) X 

— (rA,a;) 
x^ 



representing the dimensionless ground-state-atom vdW energy in terms of the trans- 
verse and longitudinal rA-dependent (local) photonic DOS functions {rA,x) and 
^"(rA,x), which are given in an explicit form by Eqs. ()42() and ()43() . respectively. 
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5.1 Qualitative Analysis 

Eqs. (|68 |1 - (|7U|1 along with Eqs. (|42|1 and (|43|1 describe the ground-state-atom vdW 
energy near an infinitely long single-wall carbon nanotube in terms of the local 
(distance-dependent) photonic DOS. Eq. (|7U|) is universal in the sense that it covers 
both strong and weak atom-field coupling regimes defined as those violating and 
non-violating, respectively, the applicability domain of the conventional stationary 
perturbation theory.^ 

5.1.1 Weak Coupling Regime 

If, when the atom is faraway from the CN, the local photonic DOS .^"'"(ryi, x) is such 
small that 



e {rA,x)^l 



(71) 



and 



(r^)! <C XA, 



(72) 

then Eq. H70|) yields a well-known perturbation theory result (see, e.g., [291 l3Up of 
the form 







XA 
2tt 



2? ' ^ 



r ( \ i^i'^A^x) + ( — ) £}\rA,x) 
ro(rc) \XAj 

X + XA 



(73) 



where the first term comes from the unperturbed Hamiltonian ()24() . ()25() and the 
second one is the second order correction due to the perturbation p6|l. This can be 
equivalently rewritten in the form 



^vw{^a) 



XA 
27r 



fo(x) 



oo 

dx ■ 

x{x + xa) 



C {^A,x) ^ {rA,x) 

XA 



(74) 



from which, in view of Eqs. (|42|1 . (|43j) and basic properties of the modified cylindric 
Bessel functions (see, e.g. [SBl, fSSJ, one can immediately come to an interesting 
conclusion. Namely, if the atom is fixed outside (inside) the CN in a way that the 
weak atom-field coupling regime is realized, i.e. far enough from the CN surface, then 
the modulus of the atom-nanotube vdW energy increases (decreases) with the CN 
radius. The conclusion is physically clear as the effective atom-nanotube interaction 

^One should make a clear difference between this definition and that where the strong (weak) 
atom-field coupling regime is defined as that violating (non-violating) the applicability domain 
of the Markovian approximation I16| . Our definition here is based on the Born approximation 
associated with the strength of the coupling between the atom and the photonic reservoir, whereas 
the Markovian approximation is related to memory effects of the photonic reservoir. 
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area is larger (smaller) for larger-radius nanotubes when the atom is outside (inside) 
the CN. Important, however, is that this obvious conclusion is, strictly speaking, 
only valid in the weak atom-field coupling regime for the outside atomic position, 
while for the inside atomic position it represents a general effect of the effective 
interaction area reduction with lowering the CN surface curvature. 

In the large CN radius limit, Eq. 1)74(1 along with Eqs. (|42|) and ()43|) can be 
shown to reproduce a well-known Casimir-Polder result ^60^ for an atom near an 
infinitely conducting plane (see Appendix C). 

5.1.2 Strong Coupling Regime 

When the atom is close enough to the CN surface, the local photonic DOS ^^(r^, x) 
is large, so that one might expect the condition 



to hold true. Under this condition, Eq. ((7n|l reduces to an integral equation of the 
form 



valid in the strong atom-field coupling regime. 

Eq. ((75|) corresponds to the inverted unperturbed atomic levels described by the 
Hamiltonian (|25l) with oja given by Eqs. (|53|) and H44() (see Fig. [7|as an example). 
The levels are degenerated when the left-hand side of Eq. ()75() is exactly 0.5. A note- 
worthy point is that if one used the weak coupling approximation of Eq. (|7U() for 
the atom close to the CN surface where Eq. ((7T|) is not satisfied and Eq. ((7^ holds 
true instead, then the result would be divergent at the lower limit of integration. To 
avoid this fact, one has to use either the strong coupling approximation (|7H|) or the 
exact equation (f7n|l to calculate the vdW energy of the atom close to the nanotube 
surface. 

5.2 Numerical Results and Discussion 

Using Eqs. (|nS)) -(|7n |) along with and we have simulated the total ground 
state energy of the whole system "two- level atom + CN- modified vacuum electro- 
magnetic field" and the vdW energy of the two-level atom nearby metallic and semi- 
conducting carbon nanotubes of different radii. The dimensionless free-space spon- 
taneous decay rate in Eq. (jTO)) was approximated by the expression ro(a;) ~ 137~^a; 
valid for atomic systems with Coulomb interaction |49j . 




(75) 




(76) 



27 



Figure |7| shows the dimensionless ground state energy level of the whole sys- 
tem [given by the corresponding eigenvalue of the total Hamiltonian H23|) - H26() and 
represented by Eqs. ()68() - (|7U() ] and two dimensionless energy levels of the unper- 
turbed Hamiltonian PH) . as functions of the atomic position outside the (9,0) 
CN. As the atom approaches the CN surface, its unperturbed levels come together, 
then get degenerated and even inverted at a very small atom-surface distance. In 
so doing, the weak coupling approximation for the ground state energy [given by 
Eqs. H68() . H69() . ((71)) and (02), (031)] diverges near the surface, whereas the strong cou- 
pling approximation [Eqs. ((HH|) . (|(i9|). (|7H|) and (|1^ . (flH]) ] yields a finite result. The 
exact solution reproduces the weak coupling approximation at large and the strong 
coupling approximation at short atom-surface distances, respectively. Close to the 
nanotube surface, the exact solution is seen to be a little bit lower than that given 
by the strong coupling approximation ((76)) . This is because the condition (|75|) turns 
into the inequality in this region with the left-hand side larger than unity that is 
not taken into account by the equation ((76|) . 

The degeneracy of the unperturbed atomic levels is the consequence of the trans- 
verse local photonic DOS increase [see Eq. H53() ] as the atom approaches the CN sur- 
face. A typical example is shown in Fig. |Sl(a) for the atom outside the (9,0) CN. The 
transverse DOS function is seen to increase with decreasing the atom-CN-surface 
distance, representing the atom-field coupling strength enhancement with surface 
photonic modes of the nanotube as the atom approaches the nanotube surface. The 
longitudinal local photonic DOS behaves in a similar way as is seen in Fig. |H1 (a).^ 
The vertical dashed lines in Fig. |H1 (a) indicate the 'bare' atomic transition frequen- 
cies XA for which the vdW energies shown in Fig. |H1 (b) were calculated. Some of 
them are the positions of peaks, the others are the positions of dips of the local 
photonic DOS. They are typical for some rear-earth ions such as Europium, or for 
heavy hydrogen-like atoms such as Caesium (supposed to be non-ionized near the 
CN). For example ^18j, for a well-known two-level dipole transition ^i<2 <-> ^Dq of 
Europium at A = 615 nm one has xa ~ 0.37, whereas for Caesium one obtains from 
its first ionization potential [SJ the estimate xa ~ 3.89 eV x 3/4 x (270)^"'^ ~ 0.5 
(the factor 3 /4 comes from the Lyman series of Hydrogen) , or less for highly excited 
Rydberg states. 

Figure IHl(b) shows the vdW energies of the atom outside the (9,0) CN for three 
different 'bare' atomic transition frequencies indicated by vertical dashed lines in 
Fig. El (a). Here, exact numerical solutions of Eqs. (|68|) -(|70 |) . ((l2|) . ((i3|) are com- 
pared with those obtained within the weak coupling approximation (|74() . (|42|) . (|43() 

''For the small atom-CN-surface distances illustrated here, the electromagnetic field retardation 
effects play no role and the transverse and longitudinal DOS functions are actually exactly the same. 
Indeed, taking the limit c — > 00 allows one to set v — h in Eqs. I)f 14^ and H115|l in Appendix B, 
thereby making the transverse and longitudinal DOS functions identical in the non-retarded limit. 
Both of these functions are proportional to c"^, but this large factor is cancelled out by the factor 
in the denominator of the dimensionless vacuum spontaneous decay rate that multiplies the 
photonic DOS functions in all the expressions. 
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Figure 8: Longitudinal local photonic DOS's (a) and ground-state vdW energies (b) 
as functions of the atomic position for the two-level atom outside the (9,0) CN. The 
('bare') atomic transition frequencies are indicated by dashed lines in Fig. |Hl(a). 



29 




Figure 9: Longitudinal local photonic DOS's (a) and ground-state vdW energies (b) 
as functions of the atomic position for the two-level atom nearby the (10,10) CN. 
The ('bare') atomic transition frequencies are indicated by dashed lines in Fig.|^(a). 
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Figure 10: Ground-state van der Waals energy of the two-level atom positioned at a 
fixed atom-surface distance [equal to the radius of the (9,0) CN] inside and outside 
"zigzag" (m,0) CNs of different radii, as a function of the nanotube index m. 



and within the strong coupling approximation 1)76(1 . I|4'2j) . 1)43(1 . At small atom-CN- 
surface distances, the exact solutions are seen to be fairly well reproduced by those 
obtained in the strong coupling approximation, clearly indicating the strong atom- 
field coupling regime in a close vicinity of the nanotube surface. The deviation from 
the strong coupling approximation increases with transition frequency xa, that is 
easily explicable since the degeneracy condition ()75(1 of the unperturbed atomic lev- 
els is more difficult to reach for larger inter-level separations. As the atom moves 
away from the CN surface, the exact solutions deviate from the strong coupling 
solutions and approach those given by the weak coupling approximation, indicating 
the reduction of the atom-field coupling strength with raising the atom-surface dis- 
tance. The weak-coupling solutions are seen to be divergent close to the nanotube 
surface as it should be because of the degeneracy of the unperturbed atomic levels 
in this region. 

Figure|ni(a) shows the longitudinal local photonic DOS's for the atom inside the 
(10,10) CN at three representative distances from the nanotube wall. As the atom 
approaches the wall, the longitudinal DOS function increases in a similar way as that 
Fig.|Sl(a) for the outside atomic position. The vertical dashed lines indicate the 'bare' 
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atomic transition frequencies for which the atomic vdW energies in Fig. |^ (b) were 
calculated. Figure|Hl(b) represents the exact solutions of Eq. ((70]) for atom inside and 
outside the (10,10) CN. When the atom is inside, the vdW energies are seen to be 
in general lower than those for the atom outside the CN. This may be attributed to 
the fact that, due to the nanotube curvature, the effective interaction area between 
the atom and the CN surface is larger when the atom is inside rather than when it 
is outside the CN. This, in turn, indicates that encapsulation of doped atoms into 
the nanotube is energetically more favorable than their outside adsorption by the 
nanotube surface - the effect observed experimentally in Ref. [SI ■ 

Comparing Fig.|Sl(b) with Fig. 1^1 (b) for the outside atomic position, one can see 
that, if the atom-CN-surface distance is so large that the weak-coupling approxima- 
tion is good, then for the same distance and approximately the same atomic transi- 
tion frequency xa ~ 0.3 the vdW energy near the (10,10) nanotube {Ren = 6.78 A) 
is lower than that near the (9,0) nanotube [Ren = 3.52 A). This is in agreement 
with the qualitative predictions formulated above in analyzing the weak coupling 
regime. In the present case, the effective atom-CN-surface interaction area is larger 
for the atom outside the larger-radius (10,10) nanotube, thereby explaining the 
effect. Obviously, one should expect an opposite effect for the inside atomic posi- 
tion. Experimentally, the heat of absorbtion of Krypton, to take an example, has 
been shown to be larger for graphite (0.17 eV/atom) than for multiwall nanotubes 
(0.12 eV/atom) j62| . in agreement with our theoretical predictions. 

A similar tendency is demonstrated in Fig llOl Here, the vdW energies calculated 
from Eqs. (|E5|) - H70|) . are represented for the atom positioned at a fixed dis- 

tance from the nanotube wall [chosen to be equal to the radius of the (9,0) CN] inside 
and outside the "zigzag" (m, 0) CNs of increasing radius Ren = m (a\/3/27r). When 
the atom is inside, the energy absolute value goes down with Ren, representing a 
general tendency of the effective interaction area reduction with lowering the CN 
surface curvature. When the atom is outside, the effect would be the opposite one if 
the atom were weakly coupled to the field, thus indicating that the atom-field cou- 
pling regime is not actually weak for the distance chosen. In general, the inside and 
outside vdW energies approach each other with Ren as it should be since the two 
atomic positions become equivalent in the plane limit where Ren ^oo. Interesting 
is also the fact that the vdW energies of the atom outside the metallic nanotubes 
(m is divisible by 3) are on average a little bit smaller than those for the atom 
outside semiconducting ones. This is certainly the property related to the difference 
in conductivities of metallic and semiconducting "zigzag" nanotubes. The property 
can be understood qualitatively as follows. In view of the absence of the longitudinal 
depolarization [lU], the CN dielectric tensor zz-component ezz(R,^) is related to 
the longitudinal polarizability azz{^) of the CN per unit length via the equation 

e^2(R, w) = 1 + iirpTlrazzii^) 
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with It being the length of the tubule. This, via the Drude relationship Q, yields 



from which it immediately follows that larger real conductivities inherent to metallic 
CNs lead to larger imaginary polarizabilities corresponding to stronger optical ab- 
sorbtion by metallic CNs compared with semiconducting ones. It is then obvious that 
virtual long-wavelength photon exchange responsible for the long-range London-type 
dispersion vdW interaction is suppressed for the atom outside the metallic CNs due 
to their stronger absorbtion compared with the semiconducting CNs. 

6 Conclusion 

We have developed the quantum theory of near-field electrodynamical properties of 
carbon nanotubes and investigated spontaneous decay dynamics of excited states 
and van der Waals attraction of the ground state of an atomic system (an atom 
or a molecule) close to a single-wall nanotube surface. In describing the atom-field 
interaction, we followed the electromagnetic field quantization scheme developed for 
dispersing and absorbing media in Refs. [12i ll.Sj . This quantization formalism was 
adapted by us for a particular case of an atom near an infinitely long single- wall CN. 
We derived the simplified secondly quantized Hamiltonian representing the " atom- 
nanotube" coupled system in terms of only two standard approximations. They 
are the electric dipole approximation and the atomic two-level approximation. The 
(commonly used |25j ) rotating wave approximation was not applied and the diamag- 
netic term of the atom-field interaction was not neglected. Starting with this general 
Hamiltonian, we have obtained an evolution equation for the population probabil- 
ity of the upper state and a vdW energy equation of the lower (ground) state of 
the two-level atomic system coupled with the CN modified vacuum electromagnetic 
field. The equations are represented in terms of the local photonic DOS and are 
valid for both strong and weak atom-field coupling regime. 

By solving the evolution equation of the upper state of the system numerically, 
we have demonstrated a strictly non-exponential spontaneous decay dynamics in the 
case where the atom is close enough to the CN surface. In certain cases, namely when 
the atom is close enough to the nanotube surface and the atomic transition frequency 
is in the vicinity of the resonance of the local photonic DOS, the system exhibits 
vacuum-field Rabi oscillations - a principal signature of strong atom-vacuum-field 
coupling. This is the result of strong non-Markovian memory effects arising from 
the rapid frequency variation of the photonic DOS near the nanotube. The non- 
exponential decay dynamics gives place to the exponential one if the atom moves 
away from the CN surface. Thus, the atom-vacuum-field coupling strength and the 
character of the spontaneous decay dynamics, respectively, may be controlled by 
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Figure 11: Schematic of the energy levels and the absorbtion line-shapes expected 
in the optical absorbance experiment with atomically doped CNs. In the center 
are the unperturbed atomic levels given by the Hamiltonian (|25j) . On the left and 
on the right are the levels of the coupled " atom-nanotube" system in the weak and 
strong atom-vacuum-field coupling regime, respectively. The system is excited by an 
external optical radiation. The allowed and forbidden (crossed) optical transitions 
are shown by vertical dashed arrows. See Eqs. (|53|) . (|62|l and (|44j) for Hja and O. 



changing the distance between the atom and CN surface by means of a proper 
preparation of atomically doped CN systems. 

We would like to emphasize a general character of the conclusion above and its 
applied and fundamental significance. We have shown that similar to semiconductor 
microcavities j^l^llH] and photonic band-gap materials ^lUni) carbon nanotubes 
may qualitatively change the character of the atom-electromagnetic-field interaction, 
yielding strong atom-field coupling. The study of such phenomena was started awhile 
ago in atomic physics [HHl and still attracts a lot of interest in connection with various 
quantum optics and nanophotonics applications ^2 [23 1^3 ESI ^ ^^^^ ^ quantum 
computation and quantum information processing [HI I66| 16 7| I68j . The fact that the 
carbon nanotube may control atom-electromagnetic-field coupling opens routes for 
new challenging applications of atomically doped CN systems as various sources of 
coherent light emitted by dopant atoms. 

Strong atom-vacuum-field coupling we predict will yield an additional structure 
in optical absorbance/reflectance spectra (see, e.g., |131 E]) of atomically doped 
CNs in the vicinity of the energy of an atomic transition. The illustration is shown 
in Fig. 111! Weak non-Markovity of the atom-vacuum-field interactions (that yield- 
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ing non-exponential spontaneous decay dynamics with no Rabi oscillations) will 
cause an asymmetry of an optical spectral line-shape (not shown in Fig. 111(1 sim- 
ilar to that taking place for the exciton optical absorbtion line-shape in quantum 
dots ,69, . Strong non-Markovity of the atom- vacuum- field interactions (yielding non- 
exponential spontaneous decay dynamics with fast Rabi oscillations) originates from 
strong atom-vacuum-field coupling with the upper state of the system splitted into 
two "dressed" states (see Fig. ^2 on the right). This will yield a two-component 
structure of optical absorbance/reflectance spectra similar to that observed for exci- 
tonic and intersubband electronic transitions in semiconductor quantum microcav- 
ities 051 ES]. 

The atomic state strongly coupled to the surface vacuum photonic modes of the 
nanotube is nothing but a 'quasi-lD cavity polariton' similar by its physical nature to 
quasi-2D excitonic jlH]) intersubband electronic [IHI and quasi-OD excitonic |19| I21j 
polaritons in semiconductor quantum microcavities as well as to quasi-OD excitonic 
polaritons in photonic crystal nanocavities [22] recently observed experimentally. 
Obviously, the stability of quasi-lD atomic polaritons in CNs is basically deter- 
mined by the atom-nanotube van der Waals interaction. This interaction, however, 
originates from strong atom-field coupling and, therefore, cannot be correctly de- 
scribed in terms of vacuum-QED-based (weak-coupling) vdW interaction models 
as well as in terms of those based upon the linear response theory. We have de- 
veloped a simple quantum mechanical approach to the ground-state vdW energy 
calculation of a two-level atomic system near a CN. The approach is based upon the 
perturbation theory for degenerated atomic levels, thus accounting for both weak 
and strong atom-field coupling and thereby covering the vdW interactions of the 
ground-state quasi- ID atomic polaritons in CNs. Within the framework of this ap- 
proach, the vdW energy is described by the integral equation represented in terms 
of the local photonic DOS. By solving this equation numerically, we have demon- 
strated the inapplicability of weak-coupling-based vdW interaction models in a close 
vicinity of the nanotube surface where the local photonic DOS effectively increases, 
giving rise to an atom-field coupling enhancement followed by the degeneracy of the 
unperturbed atomic levels due to the diamagnetic interaction term. 

The fundamental conclusion above is supplemented by those important for vari- 
ous applications of atomically doped CNs in modern nanotechnology. In particular, 
we have studied CN surface curvature effects on the atom-nanotube vdW interac- 
tions and have shown that an inside encapsulation of doped atoms into the nan- 
otube is energetically more favorable than their outside adsorption by the nanotube 
surface, in agreement with the experimental observations reported in Ref. 0. More- 
over, if the atom is fixed outside the CN in a way that the weak atom-field cou- 
pling regime is realized, i.e. far enough from the CN surface, then the modulus of 
the atom-nanotube vdW energy increases with the CN radius because the effective 
atom-nanotube interaction area is larger for larger-radius nanotubes in this case. For 
inside atomic position, the modulus of the atom-nanotube vdW energy decreases 
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with the CN radius, representing a general effect of the effective interaction area 
reduction with lowering the CN surface curvature. 

Finally, we here only dealt with the simplest manifestations of the strong cou- 
pling regime in the atom-electromagnetic-field interactions near carbon nanotubes. 
Similar manifestations of strong atom-field coupling may occur in many other atom- 
electromagnetic-field interaction processes in the presence of CNs, such as, e.g., 
atomic states entanglement [201, interatomic dipole-dipole interactions [ZDESl, cas- 
cade spontaneous transitions in three-level atomic systems [7H], etc. A further in- 
triguing extension of the present work could also be the study of the vdW interac- 
tions of excited atomic states where, even in the weak atom-field-coupling regime 
and in the simplest case of an atom near a planar semi-infinite medium, very inter- 
esting peculiarities (e.g., oscillatory behavior) were recently shown to exist |3Uj . 
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A Two-level approximation 

In this Appendix, we rewrite the Hamiltonian p8j) - (|22)) in terms of a two-level 
atomic model [25] . Within the framework of this model, the spectrum of the atomic 
Hamiltonian (|2n|) is approximated by the two eigenstates, upper \u) and lower \l), 
with the energies hwu and hioi, respectively. One has then 

Ha ^ huJu\u){u\ + huji\l){l\ (77) 

with a completeness relation 

\u){u\ + \l){l\= i. (78) 

Subtracting a constant term (h/2){uJu+uJi)I from the Hamiltonian (|77|) and thereby 
placing the energy zero in the middle between the two energy levels, one arrives at 
the 'bare' two-level atomic Hamiltonian of the form 

HA = ^<r, (79) 

with az = \u) {u\ — \l) {l\ and uja = oJu—oji being the 'bare' atomic transition frequency. 

In terms of such a two-level scheme, the atomic dipole moment operator d has 
the matrix elements (n|d|/) = (/|d|u) = d and (u|d|ti) = (/|d|/) = 0. Keeping this 
and the completeness relation (|78jl in mind, the interaction Hamiltonians H21|) and 
(|22l) can be rewritten as follows. 
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For the interaction H21|) . using Eqs. Q and @ and a well-known quantum 
mechanical operator equality 



{■pi)a = mi — (vi)c 



ih 



[{ri)a,HA\ , 



(80) 



where the Latin and Greek indexes enumerate particles in the atom and vector 
components, respectively, and Ha is given by Eq. ((7^1), one has 



duj — E (r^,cj)— h.c. 
^ 



E (r^, u) + h.c. 



(81) 



with a = \l){u\ and o"^ = \u){l\. Using further Eqs. (fTU]) . (|T^ and © and assuming 
the longitudinal (along the CN axis) atomic dipole orientation due to the dominant 
axial polarizability of the nanotube |401 1411 WR I43L I44| . one arrives at the secondly 
quantized interaction Hamiltonian (|26j) with the interaction matrix elements (|29|) 
and dSni). 

For the interaction (|22|) . one may proceed as follows 



H 



(2) 
AF 



EE 



5ij5af3Aa{YA)Af3{rA) 



i,j a,f3 

nii 



2fe2 



E—[{Vi)a,dp\Aa{rA)Afi{YA) 
m ■ 



(82) 



where the product of the two Kronecker-symbols in the first line is represented in 
terms of the 'coordinate-momentum' commutator as 



SijSa0 = [iPi)a, (rj)/3]. 

Using Eq. (|8U() for the momentum operator components and the completeness rela- 
tion ((75)) . one arrives at 



H 



(2) 



LOA 



d • A(r^) 



which, upon substituting Eq. (jS)) for the vector potential operator, becomes 



(83) 
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+ h.c. } . 
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The four items here can, in view of Eqs. ()1U() . (jlbf) and ©, be classified as fol- 
lows. The first item and its hermitian conjugate describe the processes with simul- 
taneous annihilation and creation of two photons (two-photon transitions). They 
may be safely neglected if the field intensity is low enough that, we believe, is the 
case for the vacuum electromagnetic field we deal with in considering the vdW inter- 

" (2) 

actions. Furthermore, keeping them in the Hamiltonian Hj^p is meaningless as the 
Hamiltonian H^p is nothing but a dipole approximation which, by its definition, 
neglects two-photon transitions. The second item and its hermitian conjugate are 
combined together by using bosonic commutation relations ((2) to give in terms of 
the bosonic field operators and the transverse dipole interaction matrix element pUf) 
the expression 



g-^Ly"dR|g^(rA,R,w)|' 
J dwdw'y (fRffR' [g^(rA,R,w)J g^(rA, R', u;')/^(R-> ^0 



2a 



f^A Jo 

where the second term is nothing but a two-photon correction to the dipole inter- 
action (|26j) . that can be easily seen by noting that az = a^ a — aa^ . This correction 
must be neglected for the reasons just discussed, so that finally one arrives at 

H^^l--^J2oJdR\g^{rA,K,io)\\ (84) 

Eq. (|84|) can now be combined with the 'bare' atomic Hamiltonian (|79|) to yield 
the 'effective' unperturbed atomic Hamiltonian H25|) with the renormalized atomic 
transition frequency given by Eq. (|28)) . 

Thus, we have proved that the Hamiltonian H18() - (|22l) . being approximated in 
terms of the two-level atomic model, is represented by the Hamiltonian (|23|) - (|2f 



B Green tensor and local photonic DOS 

(a) Green tensor of a single-wall carbon nanotube 

We start with Eq. H17|) for the Green tensor of the electromagnetic subsystem. This 
equation is a direct consequence of the equation 

A'jr 

^ (Vx Vx -e)^^E,{r,u;) = i—kl{v,u:) (85) 

obtained by substituting the magnetic field operator from Eq. (jl^f) into Eq. (|14() . 
On the other hand, under the condition of the Coulomb gauge (see, e.g., [SI|), one 
has 

E(r,Lj) = E"^(r,u;) +E"(r,Lj) (86) 
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with [compare with Eqs. © and ©] 



E {r,uj) = ikA{^,uj) , (87) 

where 



V • A(r,a;) =0, 
and 

^\r,oj) = -Vip{r,oj) . (89) 

From Eqs. (|5S)) -()89 |) . in view of Eq. (fTU]) and ((TB|) . one obtains the fohowing equation 
for the Green tensor components 

^ (Vx Vx [^C(r,R,u;) + IIC(r,R,u;)] = 5(r - R) (90) 

a=r,ip,z 

with additional constraints 

^Va^Ga,ir,R,u;) = (91) 

Q=r,ip,z 

and 

e„^7V;3"G^.(r,R,w) =0, (92) 

/3,'y=r,ip,z 

where Cap^ is the totally antisymmetric unit tensor of rank 3. Keeping Eqs. ()91() 
and (|92|) in mind, Eq. (|90|) is rewritten to give two independent equations for the 
transverse and longitudinal Green tensors of the form 

(A + A:2)^G,,(r,R,u;) = -5i{T-K), (93) 
A;2|lG,,(r,R,a;) = -4(r-R) (94) 

with the transverse and longitudinal (5-functions given by Eqs. and 1)12(1 . re- 
spectively, and -'"^"^^^^(r, R, w) defined by Eq. (|31|) . 

Electromagnetic properties of the " atom-nanotube" coupled system are, ac- 
cording to Eqs. (|H ^ - (|34|) . determined by the Green functions ■^G^^{rA,rA,uj) and 
,rA,io). To derive them, we start with -'-^^^(r, r^i, w) for which we use the 
differential representation 

^C(r, TA, oj) = V,V, + 1^ g{r, r^, w) (95) 

equivalent to Eq. (|3H) . Here, g{r, r^, w) is the Green function of the scalar Helmholtz 
equation which satisfies the radiation condition at infinity and boundary conditions 
on the CN surface. Indeed, substituting Eq. ((95|) into Eq. ((93]) and using Eq. ((9T|) . 
one straightforwardly obtains 

(A + k^)g{r,rA,to) = -6iir-rA) (96) 
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- the scalar Helmholtz equation with a transverse (5-source. For the longitudinal 
Green function llG'22(r, r^i, w), one analogously has 

^^Gzzir,rA,uj) = --^\7zVzg{r,rA,uj), (97) 
which upon substituting into Eq. (jMI) yields 

VzVzg{r,rA,uj) =6l{r-rA). (98) 
Eq. ()96() has a known solution 

5o(r, TA, ^) = -r 1 r (99) 

47r |r — r^l 

satisfying the radiation condition at infinity (see, e.g., [IH])- In our case, however, 
the functions 

^Gaz{r,rA,oj) = Q^Vc^Vz + Saz^ g{r,rA,u;) (100) 

and g{r, r^i, w), respectively, are imposed one more set of boundary conditions. They 
are the boundary conditions on the surface of the CN. Using simple relations 

E^{r,uj)=ik^Gc,z{r,rA,uj) (101) 
H^{r,u;) = -^ ^ eo^p^VpE^{r,u;), (102) 

f5,-y=r,ip,z 

valid at r 7^ for a classical electromagnetic field under the Coulomb-gauge con- 
dition [SJI, they can be derived from the classical electromagnetic field boundary 
conditions of the form 

E^\r=Rc„+0 - E^\r=R^„-0 = , (103) 



EJr=R,r^+0 - Ejr=Rcn-0 = , (104) 
Eip\r=Rc„+0 - H^\r=R^„-0 = — (^zz{Rcn, ^)Ejr=R^„, (105) 
Ez\r=R,„+0 - Hjr=R^„-0 = (106) 

(spatial dispersion neglected) obtained in Ref. 

Let VA > Ren (the atom is outside the CN) to be specific. Then, the function 
g{r,rA,uj), being the complete solution of Eq. 1)96(1 . is represented as a sum of 
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a particular solution of the inhomogeneous equation and a general solution of the 
homogeneous equation in the form 



g{r,rA,iv) = i . „ (107) 

[ >{r,uj), r < Ren 

where go{r,rA,i^) is the point radiative atomic source function defined by Eq. 1)99(1 
and (7^^-' (r, lo) are unknown nonsingular functions satisfying the homogeneous Helm- 
holtz equation and the radiation conditions at infinity. We seek them using integral 
decompositions over the modified cylindric Bessel functions Ip and Kp as follows |51j 

Jc\Bpih)Ipivr) J ^ ' 



p= — OD 



and 



1 °° r 

5o(r, rA,u;) = -— ^ V e^^'^ / Ipivr) Kp{vrA) e'^'dh , > r , (109) 



p=— oo 



where Ap{h) and Bp{h) are unknown functions to be found from the boundary condi- 
tions (Uni-dinni) inviewof Eqs. (|Tnni) . (fTTlT]) and (UHl, v = v{h,uj) = ^h? - k'^. The 
integration contour C runs along the real axis of the complex plane and envelopes 
the branch points ib/c from below and from above, respectively. The coordinate 
system has been fixed as is shown in Fig. ^ 

The boundary conditions (|TnHl) - (|TnHl) with Eqs. (fTTini) . ([TITT|) and (111^1) taken 
into account yield the following two independent equations for the scalar Green 
function ((TU7|) 



go{r,rA,i^)\r=R,„ + g^^HT^,^)\r=R,„ =9^ Hr,^ 



--Ra 



dg(+\r,u;)^ dg^-\v, 



Qj. \r—Rcn \r—Rcn 

I flr \( , I,2^ (-)( \\ 5ffo(r,rA,u;)| 

+ j9^ ^{v,u:)\r=R,^= — |r-=R,„, 

where P{uj)=^7ri azz{Rcn,^)/^- Substituting the integral decompositions ()108() and 
(|109|) into these equations, one obtains the set of two simultaneous algebraic equa- 
tions for the functions Ap{h) and Bp{h). The function Ap(h) we need (we only need 
the Green function in the region where the atom is located) is found by solving this 
set with the use of basic properties of cylindric Bessel functions (see, e.g., jKHl 1^). 
In so doing, one has 

RcnP{uj)v^I^{vRcn)Kp{vrA) 

Ap[h) - 



(27r)2[l + /5(a;) V^R,nIpivRcn)KpivR,n)] 
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These Ap{h), being substituted into Eq. (|1U8() . yield the function g^^\r,io) sought. 
The latter one, in view of Eq. H1U7() . results in the scalar electromagnetic field Green 
function of the form 

g{r,rA,uj) = go{r,rA,uj) 

(2vr)l£- Jc 1 + /?(^) V^RcnIp{vR,n)Kp{vR, - ^ ' ^ ^ 



p= 



where r^i > r > Ren- One may show in a similar way that the function g[r, r^, w) for 
r <rA< Ren is obtained from Eq. H11U|) by means of a simple symbol replacement 
Ip^Kp in the numerator of the integrand. 

Knowing (/(r, r^i, w), one can easily calculate from Eq. UlUUj) the components of 
the electromagnetic field Green tensor -'-GQ,^(r, fa, cj) and, consequently, 
^G^^(r,rA,u;) and "^^^(r, r^, cj) defined by Eqs. and respectively, which 
we actually need. 

(b) Local photonic DOS near a single-wall carbon nanotube 

Using Eqs. (PT|) and the property g{r, r^, w) = g{rA, r, u) which is obvious 

from Eqs. (|99|) and (|11()|) . it is not difficult to prove that 

^G^^{rA,rA,uj) = ■^Gzz{rA,rA,uj) , (111) 

,rA,uj) = "Gzz(rA,'rA,^) , (112) 

^GUrA,rA,u;) = 0, (113) 

whereupon, making use of the definition ^ and Eqs. (l3fil)-(BTIl. one arrives for 
rA > Ren at the r^-dependent transverse and longitudinal local photonic DOS 
functions of the form 

-± 3R^ ^ f dhl3{u:)vHl{vR,n)Kl{vTA) 

^ ^'■^'''^ 27rA;3 '^^£^Jcl + m^''RcnIp{vRcn)Kp{vRcn)' ^ ^ 

Th .^^_Rcn [ dh(3{u)h''v^l{vR,n)Kl{vTA) 

^ ^""^'"^^ 27,k^ '^^£^Jcl + (i{u;)v^RcnIp{vRcn)Kp{vR,n)' ^ ' 

Here, the sign in front of Eq. (|115j) has been changed to be positive. This reflects 
the fact that the right-hand sides of Eqs. H96() and (|U5|) are of opposite signs. As 
a consequence, partial solutions of the corresponding homogeneous equations should 
be taken to have opposite signs as well. This yields a correct (positive) sign of the 
longitudinal local photonic DOS which, along with the transverse local photonic 
DOS, must be a positively defined function. For rA < Ren, Eqs. H114() and (|115() 
should be modified by the replacement <-> Ren in the Bessel function arguments 
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in the numerators of the integrands. In dimensionless variables H44|) . after the trans- 
formation of the integration variable h = ky, these equations are rewritten to give 
Eqs. (EH) and (jH. 



(c) Proof of Eqs. dSSJ and 

Based upon the properties one can easily prove Eqs. and (|33j) for 

the dipole interaction matrix elements. The proof of Eq. H32[) is straightforward. In 
view of Eq. (|3U|) , one has 

/If? .2 J 2 f 
dR|g^(ll)(rA,R,a;)P = nhu;—-^ / dRRe a,,(R, cj) 

X ^(ll)G,,(r^,R,a;)^(ll)G,,(r^,R,u;)*. (116) 

Using further the integral relationship 

47r f ^ 
ImGa/3(r,r', w) = —k dRRe cr^^(R, u) Gaz{r,'R,uj) Gf3z{r' ,K,uj) (117) 

[which is nothing but a particular 2D case of a general integral relationship proven 
for any 3D electromagnetic field Green tensor in Ref. 13 , with Eq. @ taken into 
account] , one obtains for the right-hand side of Eq. H116|) 

whereupon using Eqs. (f^l]) and one arrives at Eq. H32() . 

The proof of Eq. (|33j) is based upon Eq. In view of this equation, the 

right-hand side of Eq. H33() takes the form 



(118) 



|dR|gW(r^,R,^)|2 = i^ro(..) 

± — 2Re /dRg^(r^,R,w)gll(r^,R,w)* 

where the last item can be further rewritten in terms of Eqs. H3U() . ()31|1 and ()117() 

to give 

SivuJAd^ _L I, 
^ 2 ^™ ^^^(rA.rA,^ , 

which is zero, according to Eq. (|113j) . Thus, one arrives at Eq. I|33|) . 
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C Relation of Eq. ( 1701 ) with the Casimir— Polder formula 



In deriving an 'infinitely conducting plane' result (the Casimir- Polder formula ||60j) 
from our theory, we will follow a general line of the work by Marvin and Toigo [73] 
who did the same within the framework of their linear response theory. 

We start with our weak-coupling-regime equation ()74() . First of all one has to 
simplify the local DOS functions {ta, x) in this equation by putting Wzz oo. 
From Eqs. (HH) and (flH)) . one has then for ^ Ren 



p= — oo 



-f^p[a^\/y^ - 1 - ieu{Rcn)x] 



where e is an infinitesimal positive constant which is necessary to correctly envelope 
the branch points of the integrand in integrating over y. 

The next step is taking a large radius limit in Eq. pi9|) . This can be done by 
using the relationship 

a^oo n=-oo ' 
b^oo 
(a — b = const) 

proved in Ref. |74j . One has 



- Im / dy\y "2 Ko[2i^xVy^-l-ie ] , (120) 



(.h^A,x) j 



where n = Ijol/hc with I = va — Ren being the atom-surface distance. The in- 
tegration variable y in this equation contains implicit frequency dependence. In- 
deed, comparing dimensionless equations (03), with (|114() . (|115() . one can see 
that y = h/k = h[c/uj) = h{ch/2'yQx). For the following it makes sense to ex- 
tract this frequency dependence by the substitution y = X/x with the dimensionless 
A = h{ch/2jQ). The result takes the form 

l^^^hrA,x) = ^ {/^(ll)(r^,x) -/^(ll)(r^,x)*} (121) 



with 
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Substituting Eq. (fnT|) into Eq. ((HI), one has 



3xA J dxToix) 
4ivr2 Wq {xa+ x)x'^ 

^ dxfo(x) 



/ (rA,x) /"(rA,x) 

x^ 



(123) 



{xA ~ x)x'^ 



f {rA,-x) + — PirA,-x) 



XA 



Here, in the second item, the change of the integration variable from x to —x has 
been made and it has been taken into account that ro(— x) = — ro(x). Noting further 
that both integrands in Eq. (|123jl do not have poles in the upper complex half plane, 
one may rotate the integration path of the first and the second integral by 7r/2 and 
by —it/2, respectively. Then, both integrals become the integrals over the positive 
imaginary axis of the complex plane and Eq. (|123|1 takes the form 



£vw\^A) 



3xa\ dzTojz) 
AiTT^ jio {xA+z)z^ 



f-^{rA,z) fH^A,z) 



XA 



(124) 



{xA- z)z'^ 



f^(rA,-zf+—fHrA,-zf 
XA 



Making further the change of the integration variable from z to iu with real non- 
negative u, one arrives at 



+ 



3xA j f° ° duTojiu) 
471"^ \ Jo ixA+iu)u'^ 

°° duTniiu 



lU 



f (rA,m) P{rA,iu) 

XA 



(125) 



/ {rA,-iu) H /"(ryi,-zu) 

XA 



/o (XA-IU)U 

Here, the functions f^^^^\rA,iu) have the following explicit form 



f-^{rA,iu) 
fHrA,iu) 



ui drj 



1 



1 



(126) 



obtained from Eq. ()122() after the change of the integration variable by the subse- 
quent substitutions x = iu and r]'^ = /u^ + 1. The functions f-^^^^\—iu)* = f^^^^\iu) 
as a consequence of a general property K^{z)* = Ky{z*) 56''. Substituting Eq. (|126() 
into Eq. (|125j) and making allowance for the fact that TQ{iu) = —iTQ{u), one further 
has 



£ijw(ryi) 



3xA f^dutoiu) 



27r2yo x\+v? J I 



drj 7] \/rf 



1 + 



\Koi2fiu7]), (127) 
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which upon the substitution rf = p becomes 



.™(r.) = I ^ I * (v^+ -^)/^„(2M«VP) . (128) 

In Eq. ()128p . the first integral over p is taken exactly by means of the relation- 
ship m 

dpy/i^Ko{a^) = J — K3/2{a)^ — (l + a)e-'^ 

Jl V O a^oo d 

and the second one - by making a large argument series expansion of the integrand 
and its subsequent termwise integration. This, upon returning back to dimensional 
variables by means of Eqs. (|69j) and using the definition (see, e.g., |l9l ITT] ) 

2 d^uJA 



for the polarizability tensor zz-component of the two-level system, finally yields 



h f°° 



u 



l + 2-l + 2{-] I 



2 2 



This is the half of the well-known Casimir-Polder result for the vdW energy of an 
atom near an infinitely conducting plane (see Refs. j6Ul I74j l. in agreement with the 
fact that our model of the atom-electromagnetic-field interactions in the presence of 
a nanotube only takes the longitudinal (along the nanotube axis, or, equivalently, 
parallel to the plane in the large nanotube radius limit) atomic dipole orientation 
into account. Another half of the Casimir-Polder vdW energy would come from the 
transverse (perpendicular to the plane) atomic dipole orientation which we have 
neglected. 
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